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1. INTRODUCTION 

Remotely sensed data have geometric characteristics and representation 
which depend on the type of the acquisition system used. To correlate such data over 
large regions with other real world representation tools like conventional maps or 
Geographic Information Systems (GIS) for verification purposes, or for further 
treatment within different data sets, a coregistration has to be performed. In addition 
to the geometric characteristics of the sensor there are two other dominating factors 
which affect the geometry: the stability of the platform and the topography.There are 
two basic approaches for a geometric correction on a pixel-by-pixel basis: (a) A 
parametric approach using the location of the airplane and inertial navigation system 
data to simulate the observation geometry and (b) a non-parametric approach using 
tie points or ground control points (Itten and Meyer, 1993). It is well known that the 
non-parametric approach is not reliable enough for the unstable flight conditions of 
airborne systems, and is not satisfying in areas with significant topography, e.g. 
mountains and hills. The present work describes a parametric preprocessing 
procedure which corrects effects of flight line and attitude variation as well as 
topographic influences and is described in more detail by Meyer (1993c). 

2. BASIS OF THE STUDY 

Test site and image data: The area "Zug-Buochserhom" is the standard test site 
of the Remote Sensing Laboratories, University of Zurich-Irchel in Central 
Switzerland. The region was covered by the AVIRIS flight #910705, run 6 of the 
NASA MAC Europe'91 campaign providing a data swath with an average nominal 
pixel size of about 18m. The first scene, Zug, represents a hilly area with highest 
elevation differences of about 600m and slopes with typical angles between 15° and 
60°. The second scene, Rigi, is an example of mountainous terrain with elevation 
differences of about 1400m and maximum slope angles up to 90°. 

AVIRIS auxiliary data: The quality assessment for the actual data set is 
described in detail in Meyer et al. (1993b). The current work uses the navigation data 
roll, pitch, and true heading (generated through the ER-2 Inertial Navigation System) 
and the roll and pitch of the AVIRIS instrument's precision gyros. 

Digital elevation model (DEM): The test area is covered by the two digital 
models (DHM-25) Zug and Rigi generated by the Swiss Federal Office of 
Topography . They have a resolution of 25m in i and j direction and of 0.10 m in 
elevation e with an average error in elevation of 2.2m±1.0m for model Zug, and 
4.4m±1.8m for Rigi. 

ADOUR conical radar tracking system: The ground-based immobile tracking 
radar system ADOUR is a dual antenna, dual frequency radar with a conical scan 
tracking system operated by the Swiss Air Force (Thomson-CSF, 1987). For the 
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current approach the three parameters latitude x, longitude y, and altitude z are 
used. The systematic error for elevation and azimuth is ± 0.2mrad and ±7m for 
distance with an update interval of 0.2 second. 

Ground reference information: A forest map was generated by scanning the 
green (forest) plate of the Swiss Topographic Map, scale 1:25,000, edition 1987 at 
50 pm with an Optronics 5040 Scanner. The average cartographic accuracy is about 
5.0 m. A shoreline map was produced by digitizing the same map in an ARC/Info 
using a digitizing tablet. The average (theoretical) accuracy is ± 8.0m. 

3. METHOD 

Figure 1 gives an overview of the core task for the new method for geocoding 
AVIRIS data. The basic goal is to reconstruct for every pixel the geometric situation at 
the time it was acquired with AVIRIS. This includes three major aspects. The first 
considers the flight line and attitude of the ER-2 aircraft, the second reconstructs the 
current observation geometry and the third treats the situation on the surface. This 
approach includes the three different coordinate systems (c,r) for the raw file, 
longitude x, latitude y, and altitude z together with roll ft), pitch <|>, and true heading y_ 
for the observation geometry and (i,j, e t j) for the DEM. 

Flight line and attitude of the ER-2 aircraft: The x,y, and z of the aircraft need to 
be known as a first step. For the 1991 European flight the information results from an 
unaided LTN90-116 INS navigation system with a position accuracy of 0.9 nmi/h 
(Perrin, 1993). These data are not accurate enough for the current approach. 
Therefore, ADOUR data are used as an alternative. The description of the attitude of 
the aircraft is based on the X from the navigation data and the to and<|> from the 
instrument data. 

Current observation geometry: The basic idea is shown in Figure 2 and 
described in more detail by Larson et al. (1994). The effort is to obtain the underlying 
surface out of the well-known location (=flight line) and the current attitude of the 
aircraft. The position vector X cr represents the location of pixel (c,r) in the aircraft 
coordinate system (x,y,z) at the instant the pixel was acquired by the instrument and 
for the ideal case where co = $ = y = 0°: 
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where c = pixel number of pixel (c,r) within line r of the raw image, maxP = 
maximum number of pixels per line (614), FOV = Field of View (in rad), and Z xy = 
altitude of ER-2 for the current x and y. The pixelwise calculation of the actual 
pointing direction includes correction of the panoramic distortion. The vector X c r is 
modified by rotations about the co, <f>, and x axis (vector X"' r , Figure 2). The 
transformed position vector X ' c " r is computed as follows: 
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Situation on surface: The topography causes a shift in the apparent pixel 
location, and affects the pixel size. The goal is now to find the intersection between 
the pixel location vector X"'j and the surface of the DEM. Within the neighborhood 
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around the transformed pixel location vector X"'j f a test vector X- • is searched for, 
which converges to X[”j . 
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where i^ ac n r = i-coordinate of the true nadir point, j^adir ~ inordinate of the true 
nadir point, j = elevation of the test point at position (i,j), and Z X/ y - altitude of 

ER-2 for the current x and y. To allow for a more precise selection of the 

corresponding surface point, the DEM oversampled to a grid size of 6m is used. To 

* 

define the intersection point on the surface, the normalized dot product of X- ■ and 

X' t "j is calculated. The vector X- j best representing X/'y is that for which the dot 

product DP has the smallest difference from 1. To represent the pixel size dependent 
on the topography, the four corner points (6m grid size) of every pixel are separately 
calculated. To prevent changes to the radiometric characteristics, the original value is 
selected by an improved extraction algorithm during the resampling to an 18m grid 
size, thereby eliminating the need to interpolate the values. 

4. DISCUSSION AND OUTLOOK 

There are currently no well-established methods of quantitatively assessing 
the success of a geocoding process. Visual inspection provides useful information, 
but cannot be used to intercompare methods. Statistical results based on residual 
calculation of single ground control points allow only a local error assessment. 

For the discussion, the Rigi scene is selected because of the more challenging 
topography. Figure 3 shows bands 13, 18, and 28 of the geocoded image overlaid by 
the scanned forested areas (green line) and digitized shoreline (blue line). The 
enlarged areas are selected dependent on their aspect and slope angle. In general, the 
results show a good correspondence between the geocoded image and the map for 
all existing topographical locations. A few locations show minor miscorrespondence. 
These problems are almost always restricted to single pixels and no general tendency 
can be recognized. The errors may result from changes in reality between the time 
the aerial photographs (which are the basis for the topographic maps) were acquired 
(1987) and the AVIRIS data acquistion, and on the fact that maps are the result of a 
generalization while AVIRIS displays every occurrence within its resolution 
characteristics. Figure 3B demonstrates another problem of the ground reference 
information. While the lower regions of the forested slope correspond perfectly, the 
AVIRIS image extends significantly the forest area on the upper limits. This 
"misregistration" is based on the problem of the determination of the forest border 
along the timber line and its representation in the forest map through symbolic point 
signatures which are suppressed through the scanning process. Figure 3D and Figure 
3E portray areas with rapidly changing slope angles. The blue shoreline proves for 
subarea (D) and the green forest line for subarea (E) the good correspondence. The 
additional verification calculates the average deviation of the geocoded image 
compared with the forest and lake map for checkpoints. Lines 1 and 2 of Table 1 
indicate the results of this test and confirm the visual validation. No systematic 
deviation was found. To double-check the performance of the new parametric 
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approach, scene Rigi was geocoded using the improved, non-parametric rubber-sheet 
approach (Itten and Meyer, 1993). Table 2 (line 3) presents the result for the improved 
rubber-sheet approach and shows the better performance of the parametric solution. 

The datasets Zug and especially Rigi need georadiometric corrections. 
Atmospheric corrections using radiative transfer models like MODTRAN-2a (Green 
et al., 1993) as well as a compensation for the slope-aspect-dependent illumination 
difference (Meyer et al. ,1993a) should complete the preprocessing of the sensor, 
system and scene related effects of the current data set. 

The whole procedure was implemented using the IDL (Interactive Data 
Language, a proprietary programming language. Research System Inc., 1993). 
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Figure 1: Overview of the data preparation and core tasks of the geocoding approach, 
where x=latitude, y=longitude, z=altitude of the airplane, and c=eolumn and r=row 
of the raw file. 
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Figure 2: Principal outline of the observation geometry to calculate the transformed 
pixel location vector where <t>=pitch, co=roll, and x=true heading, and x,y,z defining 
the coordinate axis for the ideal pixel location vector and x',y',z' for the transformed 
pixel location vector. 

Figure 3 (Slide 7): (A) shows the composite rendered on the digital elevation model 
with tags (B-E) for the location of the enlarged subareas of the horizontal (non- 
rendered) composite. (B)-(E) show zoom-up (~ factor 5) parts with different aspect 
and slope angles. 


Table 1: Result of the quantitative effort for the comparison between the geocoded images and the 
forest map with the RMS for the east-west direction (i) and north-south direction 0 of the digital 
elevation model. The check points used for the non-parametric approach for scene Rigi are a 
selection from the entire number of points used for the parametric approach. 



No. of checkpoints 

RMS for i-direction 

RMS for j-direction 

Scene Zug (parametric approach) 

186 

0.07 pixel 

0.19 pixel 

Scene Rigi (parametric approach) 

309 

0.12 pixel 

0.09 pixel 

Scene Rigi (non-parametric approach) 

249 

3.57 pixel 

1.37 pixel 



